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thermodynamical consistent CJT calculation in studying nuclear matter 
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We have attempted to apply the CJT formalism to study the nuclear matter. The thermody¬ 
namic potential is calculated in Hartree-Fock approximation in the CJT formalism. After neglecting 
the medium effects to the mesons, the numerical results are found very consistent with those ob¬ 
tained from the mean field calculation. In our calculation the thermodynamical consistency is also 
preserved. 

PACS numbers: 21.65.+f, 11.10.Wx, ll.15.Tk 


I. INTRODUCTION 

Mean field theory (MFT) of quantum liadrodynamics (QHD) is very successful in describing many nuclear phe¬ 
nomena EH, Hi. The idea of MFT, in a simple description, is that at high enough nuclear density, fluctuations in 
the meson fields could be ignored, and the meson fields could be replaced by their classical expectation values or 
mean fields jj.]. However, mean fields are insufficient for a detail understanding of short-distance nuclear physics. 
The development of reliable techniques to go beyond mean field calculation is also required Q. In a renormaliz- 
able field theory, through a Green’s function formalism, one can get relativistic Hartree approximation (RHA) by 
self-consistently summing all the tadpole graphs in the baryon self-energy [Bj, Q. It turns out that RHA consists 
MFT result and the additional vacuum fluctuation corrections. Neglecting the vacuum corrections, RHA will return 
to MFT. If one further considers the meson emission and reabsorption (“exchange”) graphs in the baryon proper 
self-energy, the approximation is referred to Hartree-Fock approximation (HF) The inclusion of the exchange 

terms only makes a small correction to MFT. In fact, after renormalization to equilibrium nuclear matter properties, 
the binding energy curves in MFT and HF approximations are almost indistinguishable @]- This match between 
MFT and RHA or HF approximations indicates that MFT is equivalent to calculate the lowest order diagrams in the 
proper nucleon self-energy in field theory. Thus one would expect to calculate higher order diagrams to get the results 
beyond mean field. In j|], through an effective action formalism, the nuclear matter energy density was explicitly 
calculated at two loop lever. It was found that the higher order contributions were enormous, which altered the 
description of the nuclear ground state qualitatively. This is mainly because that an expansion in powers of loops is 
basically an expansion in the dimensionless coupling constants which are large in QHD. The quantum corrections are 
correspondingly large. So far there are no systematic, reliable calculations of field theory including higher order loop 
contributions in the study of the nuclear matter. 

When we refer to the loop expansion, we think of that in the area of studying spontaneous symmetry breaking, the 
effective action r(</>), which is the generating functional of one particle irreducible graphs, had been systematically 
calculated by summing all the relevant Feynman graphs to a given order of the loop expansion [Tcj . The effective 
action formalism had been developed in several works [n], ll 2l nTl Ibi . A notable development in this direction was 
the generalization of the effective action for composite operators initially studied by Cornwall, Jackiw and Tomboulis 
(CJT) 0- According to this formalism the usually effective action T(()>) are generalized to depend not only on 
4>{x), a possible expectation value of the quantum field <f>(;r), but also on G{x,y), a possible expectation value of 
the time-ordered product T$(x)$(y). G{x,y) is also the propagator of the field. In this case the effective action 
T((f>,G) is the generating functional of the two particle irreducible vacuum graphs. An effective potential V{(j),G), 
which is an important theoretical tool in studying the symmetry breaking and phase transition, can be defined by 
removing an overall factor of space-time volume of the effective action. Physical solutions demand minimization of 
the effective action with respect to both (f> and G [ll|. As a result the CJT effective potential should satisfy the 
stationary requirements 


and 


dV(<t>,G) 

d(f> 


= 0 


( 1 ) 


dV^,G) 

dG 


( 2 ) 


This formalism was originally written at zero temperature. Then it had been extended at finite temperature by 
Amelino-Camelia and Pi where it was used for investigations of the effective potential of the A (j> 4 theory Lg] and 
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gauge theori es 117 1. It was also applied to study the chiral symmetry breaking in effective chiral models and QCD-like 
theories El Il9t l20j| . In studying the spontaneous symmetry breaking system, like gauge symmetry breaking in the 
electroweak system M and chiral symmetry breaking in the strong interacting system fisl ] , (f> is the order parameter 
of the transition and nonzero. While in the system without spontaneous symmetry breaking, this parameter is zero. 
Thus we find the diagram expansion of the effective potential is basically the same as that of the thermodynamic 
potential in a thermodynamic system. That is to say we can use the CJT formalism to study the thermodynamic 
system like the nuclear matter. In the CJT formalism all the fields are treated as operators and it is possible to sum 
a large class of ordinary perturbation-series diagrams to infinite order that contribute to the effective potential or the 
thermodynamic potential. The form of the propagators can be consistently determined by a variational technique, 
as in ©• This resummation scheme can be used to study non-perturbative physical effects. Substantially, the MFT 
or RHA is also a certain resummation scheme. In this sense, the CJT resummation scheme can be used to make 
calculations beyond MFT. In principle the higher order contributions can be self-consistently and clearly calculated. 
As well known, MFT is thermodynamical consistent p|. To preserve the thermodynamical consistency in the study 
of the thermodynamic system in field theory is not a trivial problem. In the beyond mean field results was 
obtained by consistently calculating the one baryon loop in the proper meson self-energy. The results indicated a 
softer equation of state and a smaller compression modulus. However, the thermodynamical consistency had to be 
achieved by including additional compensatory term which needed to be properly fixed at zero temperature. In the 
CJT formalism, the thermodynamical consistency will be automatically guaranteed by the stationary conditions for 
the effective potential. However, this can be only fulfilled if the final physical results could be self-consistently worked 
out, which could be seen in our later discussion. 

Our goal is to apply the CJT formalism to study the nuclear matter. As far as we know, no such investigations have 
been found in this direction as yet. In this paper, we will use the Walecka model (also called QHD-I model) @ and 
illustrate how the thermodynamical potential can be consistently calculated in the CJT formalism. In our calculation, 
the CJT thermodynamic potential will be calculated at two-loop level with all the loop lines treated as the full 
propagators, which is called Hartree-Fock approximation in CJT formalism. However it should not be confused to the 
usual HF approximation in the study of the nuclear matter. In the CJT expansion of the thermodynamic potential, 
Hatree-Fock approximation is the lowest order in coupling constant E3- Even in this approximation, the calculation 
is still very difficult. In this paper, as the first step, we will simplify the calculation by neglecting the medium effects to 
the mesons like in the RHA or HF calculation. The thermodynamical consistency will be achieved in our calculation. 
By a numerical study, we will reproduce the binding energy curve of the nuclear matter. The effective nucleon mass 
at finite temperature and the liquid-gas phase transitions are also discussed. These results are found very consistent 
with those of MFT. 

The organization of the present paper is as follows. In section 2, we formulate the CJT formalism in studying 
the nuclear matter through the Walecka model. The gap equations are derived at Hartree-Fock approximation. In 
section 3, we solve the gap equations and thermodynamic potential consistently by neglecting the medium effects 
to the mesons. The pressure and density can be derived accordingly. We demonstrate how the thermodynamical 
consistency is ensured by the stationary conditions in our calculation. In section 4, we give our numerical results and 
make the comparisons to the results of MFT. The last section comprises a summary and discussion. 


II. CJT FORMALISM IN NUCLEAR MATTER 

We start from the QHD-I Lagrangian which can been written as 

£ = ~ m N + - gulu^)^ + ^(d^ad^a - m 2 er 2 ) - + ^m 2 (3) 

where ■*/’> & and u> are the fields of the nucleon, sigma meson and omega meson respectively and F /llf = d^u> v — d^ui^. 
According to the Lagrangian we can write down the free propagators of the nucleon, sigma meson and omega meson 
respectively as 


Go(p)= T J —, 
y- m N 

(4) 

^o(p) 2 2 ’ 

pz _ m z 

(5) 

D^{p) = - 5 -^- 5 -. 
p z — 

(6) 
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We will use the imaginary-time formalism to compute quantities at finite temperature l22l . Our notation is 

(7 > 

where (3 is the inverse temperature, f3 = 1/T. We have u> n = 2mrT for boson and oj n = (2 n+ l)irT for fermion, where 
n = 0, ±1, ±2, ■ ■ •. A baryon chemical potential /i can be introduced by replacing po = ico n with po = iui n + p in the 
nucleon propagator. 

According to the CJT formalism 0 , the expansion of the effective potential or the thermodynamic potential in 
the nuclear matter can be written as 

Q(G,A,D)=i [ In Go(p)G~ 1 (p) + i [ [G^(p)G(p) - l] - 1 f lnA 0 (p)A -\p)- 1 - [ [A^\p)A(p) - l] 

Jp Jp ^ Jp ^ Jp 

- j\nD 0 (p)D- l (jp) j [D^[p)D{p) -l]+n 2 (G,A,I>), (8) 

where G , A and D are the full propagators of nucleon, cr meson and u> meson respectively, which are determined 
by the stationary condition n 2 (G,A 7 D) is given by all the two-particle irreducible vacuum graphs with all the 
propagators treated as the full propagators. In the CJT formalism at Hartree-Fock approximation Q 2 {G, A, D) can 
be illustrated by the graphs of FiglU where the solid lines repent G(p ), the dashed line represents A (j>) and the wavy 
line represents D lll '(p). The vertices for and are — g a and gui^ respectively. The analytic expression is 

Q q 




FIG. 1: Hartree-Fock approximation to A, D). The solid, dash and wavy lines are the nucleon propagator G, the a 

meson propagator A and the ui meson propagator D respectively, q and p are the four momenta. 


n 2 (G, A, D) = J J Tr [G(q)G(q - p)A(p)} + *£ J J Tr [rG{q)YG{q - p)D ltv {p)] . (9) 

From the stationary condition © which demands that f1(G, A, D ) be stationary against variations of G, A and D 
respectively, we will have the following gap equations 

G(q)- 1 = G 0 (qy i +gl [ G(q - p) A(p) + g 2 u f yG(q-p)yD^(p), (10) 

J P J p 

A(p)- 1 = A 0 (p)- 1 -gl / Tr [G(q)G(q - p)], (11) 

Jq 

D^p)- 1 = Do^ip)- 1 - gl f Tr[^G(q)^G{q-p)}. (12) 

Jq 


The above equations can be also represented pictorially in Fig© These integration equations are nonlinear coupled 
and momentum dependent. Needless to say they are very difficult for computation. The certain approximations 
should be adopted. 


III. SOLVING GAP EQUATIONS AND THERMODYNAMIC POTENTIAL 

In this section we will solve the gap equations through certain approximations. The thermodynamic potential will 
be consistently determined. Then we will derive the pressure and the net baryon density and demonstrate how the 
thermodynamical consistency is achieved. 
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G-\q) = G 0 \q) + 

A _ 1 (p) = A 0 \p) + 



+ 



q-p 


q-p 


D 


\p) = A 1 ip) + 



Q-P 


FIG. 2: The gap equations satisfied by the nucleon, the a meson and the ui meson propagators at Hartree-Fock approximation. 
q and p are the four momenta. 


To solve the equations, we need first to decouple the equations. As in usual Hartree or Hartree-Fock approximation 
in studying the nuclear matter, the meson propagators are treated as the free propagators In this paper, as the 
first step, we will simplify the equations by replacing the meson propagators by the free ones. This also means we 
have neglected the medium effects to the mesons. We expect that under this approximation it could yield the results 
which are consistent with those of the MFT. The approximation means 

A (p) = A 0 (p), D^ip) = D 0iliV (p). (13) 

As a result the thermodynamic potential is reduced to 


H(G,A 0 ,C 0 )= i [ lnG 0 (p)G- 1 (p)+i [ [G^\p)G{p)- l] 

J p J p 

+ tJ 2 L JJ q Tr [G{q)G{q ~ P)A ° {P)] + ff p f q Tr ^ G ^ G ^ ~ P) D °MP)} ■ ( 14 ) 

The gap equations are much simplified accordingly. Only the nucleon gap equation m is left and needs to be solved. 
To proceed we take the following ansatz of the full nucleon propagator 


G(q) = 


(15) 


'y ^ - rn N + E ’ 

where E is the proper nucleon self-energy which will be determined by the equation m . It can be generally written 
as 


E(q) = E s (q) - 7 °E 0 (g) + 7 -qE„(g). 
Thus we can define an effective nucleon mass 


(16) 


M N (q) = m N - E 8 (q), 


(17) 


which is momentum dependent. As the nuclear matter is a uniform system at rest, the E„ term in equation d) is 
usually neglected as in RHA calculation fjj). In the usual HF approximation in the study of the nuclear matter, the 
contribution of the E„ term to the final physical result is also found to be very small @. In our discussion this term 
will be neglected. Substituting equations m and m into equation m, after some calculations on separating the 
terms which are associated with and without 7 matrix, and then by comparing the two sides of the equation, we can 
obtain the following equations for the components of E 


A (q) = igl [ 

J V 


A (9) = igl 


—M N (q ) 


[(95 - po ) 2 - E i- P ] \p 2 o - a(p) 2 ; 

_go - Po _ 

[(go - po ) 2 - El-pi [Po - A(p) 2 ] 


+ A 


4Mjv(g) 


■Wu. 


[(go - Po] 2 - £q-p] bo - A(p) 2 ] ’ 

_2gp* - 2p 0 _ 

[(go - Po) 2 - El_ p \ bo - A(p) 2 ] ’ 


(18) 

(19) 
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where = qo — Eg (q), E <7 / u ( p) = p 2 + m?j u and £? q _ p = ^/(q — p) 2 + Mjv(g) 2 . After we perform the sums over 
the Matsubara frequencies, we obtain 


S s{q ) = -gl j 

So(g) = glj 


d 3 p M N (q) 
(2tt) 3 2£ CT (p)£’ 1 
d 3 p 1 


r-[-^o- - -®<t] + 9u f 

q-p J 

r- [C ct -D ct ]+^ f 

q-p J 


d 3 p 2M N (q) 

(2^) 3 ^(p)A q _ p 

d 3 p 1 


[A w — _B W ] 


(2tt) 3 2£ a (p)£ q _ p L aj (2 tt) 3 £ w (p)£7 q _ p 

where the expressions of Ai, Bi, Ci and Di (i = <j,lu) are given in the following, 


[C u - D u \, 


A, = 


a = 


in which £) = -^/p 2 + : 


(do A q _ p ) fi-\-Ej 

(q* 0 -Eg-p) 2 -E 2 ■ 
<7g(<?g - ^q-p)«i - Efm - EiEq-ph+ 

(q*o - £q-p) 2 - Ef 
i 2 and 


Bi = 


D, = 


1 


^(^g A Eq—p) T fl— Ej 

(q* 0 + Eg- p ) 2 - Ef ’ 
gg(gg + Eq-p)m - Efm - EjEq- p h- 
(q* 0 +Eq_ p y-E? 


1 


Hi = 


,/SEi _ 1 ’ 


n± = 




( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


where g* = g — Eg(< 7 ) and E 4S the baryon chemical potential. In obtaining equations (1201) and EJ, we have 
encountered the divergent terms which are independent of the distribution functions rii and h±. In principle these 
divergent parts could be properly renormalized. However, in this rudimentary work, we simply neglect the divergent 
parts which are not explicit temperature dependent. This treatment is not strange in the literature [f 8l l23l |. 

Now the equations are finite but momentum dependent, which are still intractable. We note that in study ing the 
nuclear matter, the effective mass is usually defined as the pole of the full propagator in the limit q —> 0 |2j,|25|]. In 
light of this definition, the effective nucleon mass here can be defined by 


qf - q 2 ^(kol = m n, |q| = 0) 


|q|=0 


= 0. 


(25) 


This means in equations EJ and EJ we will take |q| = 0 and set |(?g| = Mn- The treatment is different from the 
usual pole approximation for that we will take |q| = 0 before the angular integration performed. This will greatly 
simplify the calculation. Furthermore, as |gj| = Mn , q q can be set to either positive or negative value. In order 
to ensure that at g = 0 the baryon density keeps zero which can be realized in our later discussion, here we will 
take = Mn in Ai and Bi, while take q g = —Mn in Ci and Di. After these approximations and treatments, the 
equations EJ and EJ will be further simplified to 


S s = ~9o 


d 3 p M n 2(M n — E)n a — E a (n + + fi-) 2 

— + g u 


d 3 p 2 M n 2(M n - E)n w - E ltJ (n + + «-) 


(27r) 3 2E a E 


So =gl 


(M n - Ef - El 


(2t r) 3 E U E 


(Mn - E) 2 - E 2 


d 3 p 


-(*4 “«-) 


(2tt) 3 2 (M n - E) 2 - 2E 2 a 


d 3 p — (n + — nJ) 

(2t r) 3 (M n -E) 2 -E 2 J 


(26) 


(27) 


where E = \/p 2 + Mjy. Now E s , Eg and Mn become momentum independent and can be numerically calculated. 

In the next we will demonstrate how the thermodynamic potential, the pressure and net baryon density can be 
derived in a thermodynamical consistent way. 

Considering equation EJ- El and El we can have 


Go(q)~ 1 G(q) - 1 = iX(q)G(q), 


Z(q) = igl f G(q-p)A 0 (p)+igl j ^pG(q - p)^D^(p). 

J p J p 

Substituting the above equations into equation EJ- the thermodynamic potential will be reduced to 

n = ij In Go(p)G(p)~ 1 -\\et [G(«z)E(g)] . 


(28) 


(29) 


(30) 
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The pressure can be obtained by the thermodynamic relation 


p = -n. 


( 31 ) 


For the consistency of the calculation under our approximation, £ (q) in equation (13011 will be also determined at 
|q| = 0 and |qg| = Mjv, which means £ becomes independent of the momentum in equation (13011 . After the frequency 
sums in equation (®) and considering equation m we can obtain the pressure as 


P(J*,T) 


4 

3 


f d 3 p p 2 

1 1 

— 2Mat£ s 

f d 3 p 1 

1 1 

/ (27 r ) 3 E 

e 0 + i ' ePiE+v*) 4 . i 

(2tt) 3 E 

e P(E-/J.*) 4- 1 1 e 0(E+v*) 4- 1 


+2Eq 


[ d 3 p 

1 

1 

/ (2tt)3 

e 0(E-li*) 4- 1 

e 0(E+fi*) _|_ ^ 


(32) 


where E s and Eo are determined by solving equations JS2> and m- In (131211 we have also neglected the divergent 
terms which are independent of the distribution functions. 

The thermodynamical consistency requires that at fixed T and p the pressure will be maximized with respect to 
E s or Eo, when E s and Eo are independent variables. This requirement will be satisfied by the stationary condition 
in the CJT formalism. To make it clear, we assume that P is still in its functional form of G , then from equation 
which indicates G is a function of E, and considering the stationary condition m we will have 


dP(G (£)) _ dP(G( S)) dG( E) _ n 

d£ dG dT, ^ ’ 

This equation ensures the thermodynamical consistency. So when the pressure or the thermodynamic potential 
is self-consistently determined by solving the gap equations, which are derived by the stationary conditions, the 
thermodynamical consistency will be automatically achieved. However, one should notice that the pressure in the 
form of equation lETl can not be maximized with respect to £ s or Eo, because in obtaining equation l!321) . the gap 
equations have been already substituted into it, which makes that E s and Eo are not independent variables in equation 
ll?T21> . 

Next we will derive the net baryon density of the system. The density will be determined by the thermodynamic 
relation 



(34) 


From a general expression of the pressure, which is a function of the independent variables T, p and E, we can write 
down the partial derivative in the following expression 


dP(E,p,T) 


dP&,p,T) 

ap(E,/i,r) 

5E 

dp 

T 

dp 

E,T ^ 

H,T 


(35) 


According to equation (ESI the second term on the r.h.s of equation (13511 becomes to zero, which shows a fulfillment 
of the thermodynamical consistency. The partial derivative is reduced to 


dP{£,p,T) 


dp 


dP&,p,T) 


dp 


(36) 


E,T 


Then from equation gd and considering equation ED, the derivative can be formally evaluated as 


dP(Z,p,T) 


dp 


,d_ 

V 


In [G 0 (p)G(pY 


Tr 


dG(q) 

dp 


S(5) 


- / Tr 


f dG(q) 
l dp 


ial [ G{q-p)A 0 (p)+igl j - P)^ D ^(p) 

J V j V 


(37) 


From equation EH> we find that the second term and third term on r.h.s of equation ED will cancel each other, thus 
the final result is 


ap(E,/i,T) 


dp 




d_ 

dp , 


In [G 0 (p)G(p)- 1 ] 


(38) 
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Then the net baryon density is given as 


P = 4 


d 3 p 
(27r) 3 


1 

e /3(B-/i*) X 


1 

e /3(iS+7t*) _|_ ^ 


(39) 


This is the standard form of net density for quasi-particles, which also indicates that the thermodynamic functions 
could be calculated thermodynamical consistently. 

Furthermore the energy density of the system can be derived by 

dP 

e=~P + pp + T—. (40) 

In the following we can study the thermodynamics of the nuclear matter. 


IV. NUMERICAL RESULTS AND COMPARISON TO MFT 


The coupling constants g G and g u will be refitted by reproducing the correct saturation properties of the nuclear 
matter at zero temperature. The nucleon mass, sigma meson mass and omega meson mass are taken as tojv = 939A leV, 
nrifj = 550 MeV and = 783 MeV. By setting T = 0, from equation E3, E! and El we can reproduce the correct 

binding energy curve of the nuclear matter which means at a saturation density of 0.16 nucleons per fm 3 it has a 
binding energy of 16 MeV per nucleon. The curve of energy per nucleon versus density is shown in FigEl The coupling 
constants thus are fixed at g 3 = 155.6 and g 3 = 521.5. These values look a little greater than those of the MFT 
which are g% = 91.6 and g% = 136.2 J6j. The binding energy curve of MFT is also plotted in FigEl in dashed line. 



FIG. 3: The average energy per nucleon minus the nucleon mass tun as functions of baryon density at zero temperature in the 
CJT formalism (solid line) and MFT (dashed line). 



e (GeV/fm 3 ) 


FIG. 4: The nuclear matter equations of state at zero temperature in the CJT formalism (solid line) and MFT (dashed line). 









We can see that the result of CJT is very close to that of MFT. At low density the energy curve of CJT is almost 
overlap with the energy curve of MFT; at high density the curve of CJT rises a little slower than that of MFT. The 
compressibility of nuclear matter at saturation density calculated in the CJT formalism turns out to be 

K = pl d2 \ £ { p) = 485 MeV, (41) 

dpp 

where pp is the Fermi momentum. It is smaller compared to that of the MFT which is K = 540MeU. The equation of 
state P vs. e for nuclear matter is shown in Fig 0| The CJT and MFT results are almost overlap. Note the approach 
from below to the casual limit e = P (where c SO und = cu g ht ) at high density. 



FIG. 5: The effective nucleon mass as functions of temperature at zero density in the CJT formalism (solid line) and MFT 
(dashed line). 



FIG. 6: The isotherms of P vs. p at T = lOMeF and T = 15 MeV in the CJT formalism (solid line) and MFT (dashed line). 

At finite temperature, from equation 0 and considering 1251) . the effective nucleon mass as a function of tem¬ 
perature at zero density can be plotted in Fig|5] (solid line). When compared to that of MFT (dashed line), the 
effective mass from the CJT calculation drops a little more slowly at high temperature. The well-known liquid-gas 
phase transition of the nuclear matter at low temperature also exists by the CJT calculation. The isotherms of P 
vs. p, which show the first order transitions, are plotted in FigEI The CJT results are again found very close to the 
MFT results. The critical temperature in the CJT calculation is about 19 MeV which is almost the same as that in 
the MFT. 


V. SUMMARY 

In this paper we have made an attempt to use the CJT formalism to study the nuclear matter. As the first step, 
we have neglected the medium effects to the mesons and obtained the results which are found very consistent with 
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those from MFT. We have also demonstrated how the thermodynamical consistency has been achieved in the CJT 
formalism in studying the nuclear matter. In our discussion one can also see that the beyond mean field calculation of 
the nuclear matter can be carried out in the CJT formalism, at least by including the medium effects to the mesons. 
However it is obvious that the calculation will be quite involved. This topic will be studied in our further work. 
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